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Abstract 

Scattering of sub-GeV dark matter (DM) particles with hydrogen atoms is studied in this pa¬ 
per. The interactions of DM with electrons and nucleons are both included and formulated in a 
general framework based on nonrelativistic effective field theory. On the assumption of same dark 
matter coupling strengths, it is found that DM-electron interactions dominate the inelastic atomic 
transitions to discrete excited states and ionization continuum around the threshold regions, and 
DM-nucleon interactions become more important with increasing energy and dominate in elastic 
scattering. The conclusion should apply, qualitatively, to practical detector species so that electronic 
and nuclear contributions in DM scattering processes can be disentangled, while issues including 
binding effects and recoil mechanism in many-body systems call for more detailed calculations. 
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I. INTRODUCTION 


The existence of dark matter (DM) has been concretely established based on gravitational 
evidences from various scales of the universe. However, its composition and non-gravitational 
interactions, if any, are still unknown. Among all DM candidates, the weakly-interacting 
massive particles (WIMPs) receive most attention mainly because of the so-called “WIMP 
miracle”: with DM particle mass ~ 10 GeV-1 TeV and self-annihilation cross section 
similar to the one of the weak interaction, WIMPs can explain the relic DM abundance of 
the universe. By looking for nuclear recoil triggered by WIMP scattering in detectors, null 
results of direct searches have ruled out quite some WIMP parameter space (for current 
status, see, e.g., [T] and references therein), and future multi-ton experiments will further 
improve the limit (or make discovery). 

On the other hand, the dark sector consisting of light DM (LDM) particles with masses 
below 10 GeV also attracts interests recently (for overview, see, e.g., |2] and references 
therein). First of all, there is no reason to exclude such possibilities a priori, and in fact, many 
well-motivated models predict their existence: for example, WIMPless and asymmetric 
DM with varying between MeV to GeV scales, several proposals at MeV ranges [MO], 
and keV ranges including bosonic super-WIMP im , axinos HMD, and sterile neutrinos 
[TU] . Furthermore, in order to accommodate low-energy anomalies in astrophysical sources 
such as the 511 keV [20l El] (reviewed in [22]) and 3.5 keV |23l El] emission lines, and the 
GeV 7 -ray excesses in the Galactic Genter [251 EE], LDM are often proposed as possible 
answers. 

As most of these LDM candidates can not produce observable nuclear recoils in present- 
day (and near-future) detectors, direct searches for them have to rely on electron recoils. 
Not like WIMP searches, they are used to constrain DM-electron interactions mm- 
Alternative searches through indirect signals and colliders are also reported jEEHS]. 

DM-electron interactions are interesting from another aspect beyond LDM: To reconcile 
the tension between the light WIMP signal (m^ ~ 10 GeV) reported by the DAMA/LIBRA 
experiment [32ll3l] and stringent constraints set by other direct experiments on WIMP- 
nucleon interactions, models of DM being lepto-philic but not hadro-philic provide viable 
solutions (see, e.g., [35H3Z]). At even higher masses, the annihilation of WIMPs are among 
the favored explanations of the positron excesses observed by PAMELA |38] l39] . Fermi- 
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LAT go], and AMS-02 gH |12| experiments. 

If DM interacts both with electrons and nncleons, an important qnestion natnrally arises: 
In a specific DM scattering process which a direct search detector is build to look for, what 
are the contributions from the electronic and nuclear degrees of freedom? Even though 
the current common practice in constraining DM interactions is one-type-at-a-time, it is 
necessary to keep in mind that events measured by a detector are a sum from all possible 
sources. Furthermore, it is desirable from the experimental point of view to determine which 
process and kinematic region would be best to constrain a certain type of DM interactions 
with electrons or nucleons. For this purpose, one has to rely on theoretical analysis. 

In this article, we attempt to address the above questions using the simplest atom; hy¬ 
drogen - where most calculations can be carried out analytically - and study its scattering 
with nonrelativistic LDM particles of a MeV-GeV mass range. As the associated energy 
and momentum transfer of such scattering processes overlap with typical atomic scales, one 
expects that atomic physics plays an important role and issues like binding effects and elec¬ 
tron/nuclear recoils require detailed study. The simplistic setup of hydrogen should therefore 
provide useful qualitative understanding which applies to more intricate systems. 

The paper is organized as follows. In Section [IT] the general formalism is developed for 
scattering channels including elastic, discrete excitation, and ionization. The DM-electron 
and DM-nucleon interactions are formulated in a general way based on effective field theory 
(EFT). In addition to commonly-discussed spin-independent and spin-dependent contact 
interactions, which are the leading-order terms in the EFT expansion, we also include the 
possibility of long-ranged DM interactions (for example, through kinetic mixing of dark and 
normal photons) and a few next-to-leading-order terms. In Section]^ we present and discuss 
our results, and infer from these generic features that will apply to other atoms including 
those practically being used in mainstream DM detectors. Finally, a summary is given in 
Section HVl 


II. FORMALISM 

Direct searches of DM look for signals as results of DM scattering off normal matter. As 
the nature of DM and its non-gravitational interactions with normal matter are still un¬ 
known, instead of considering specific, well-motivated models, we adopt a general approach 
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based on effective field theory (EFT). A nonrelativistic (NR) EFT that accommodates scalar, 
fermionic, and vector DM particles with velocity ij^and their interactions with protons 
(p) and neutrons (n) via intermediate scalar and vector bosons is formulated in Ref. |43| . 
and fully worked out to next-to-next-to-leading order in Ref. H. In this work, we further 
take electrons into account, and focus on the kinematic region where electrons also behave 
like NR particles. 

At leading order (TO), the effective interaction takes the form 


f=e,p,n 

+d'/’^M*x)(f*f)+d'r-^^(x*S:,X) ■ (PS,/)} , (1) 

where y and / denote the NR DM and fermion helds, respectively; are their spin 
operators (scalar DM particles have null S^)] the magnitude of the DM 3-momentum transfer 
g = |g| depends on the DM energy transfer T and its scattering angle 9. Note that we use the 
same nomenclature as in m for the low-energy constants (LECs) cY that characterize the 

/ f\ 

types of the y-/ contact interactions. Correspondingly the LECs d) '^’s are used to describe 

potential 17(l)-like, long-ranged y-/ interactions that are results of, e.g., mixing of dark and 

, ' 

ordinary photons via where the and refer to the held tensors of dark and 

ordinary photon, respectively, and e the mixing angle. These LECs corresponds to the ones 
of j43j by C\ —^ C 4 —^ h 2 , —y and d 4 —^ l 2 - 

To simplify the presentation of the full scattering formula for an unpolarized DM scat¬ 
tering off a hydrogen atom, we start with the simplihed case where only one of the LECs is 
assumed to be nonzero. (Note. This is the conventional practice in DM searches). Following 
the standard scattering theory (for more details, see, e.g.. Ref. |1S]), the differential cross 
section in the laboratory frame for DM being scattered by the term alone in Eq. ([^ into 
the hnal 3-momentum /c 2 with an inhnitesinial phase volume d^k 2 is expressed by 



— Y6{T - Ecu - {Ef 

F I 


El)) 


d?k2 

{2nY 


( 2 ) 


^ We adopt the natural units c = 1 and h = l. 
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The reduced mass /x = mefrip/{irie + ?Up) with rrie^p) being the mass of electron (proton); for 
later use, the mass of hydrogen is designated Mh = M — B with M = rrie + rup and B the 
binding energy. The initial state |/) denotes the hydrogen atom at the ground state, i.e., 
the spatial part |/)spatiai ~ Dirac delta function imposes the energy conservation 

that the energy deposited by DM equals to the recoil energy of the atomic center of mass, 
Eqm, plus the internal excitation energy Ep — Ej of the atom. 

Depending on the nature of the final state {E\, the scattering processes are classified as 

1. elastic scattering (el): spatial {E\ = (ls| ; Ecu = and Ep - Ej = 0, 

2. discrete excitation (ex), spatial (-^1 = (nlmi\ with {n,l,mi) ^ (1,0,0); Eqm = q‘^/{2MH) 
and Ep- El = E^i - Eu, 


3. ionization (ion): spatial (-^1 = ipA with pr denoting the relative momentum in the CM 
frame; Eqm = q^/{2M) and Ep — Ej = B +pl/{2p). 


The symbol Yli uieans an average over the initial magnetic (and spin, when spin operators 
are involved) states; YIf uieans a sum over all the final magnetic and spin (also spin, 
too) states for elastic scattering and discrete excitation, while for ionization, the sum over 
magnetic states is replaced by J d?pr/{2 ti)'^. 

The analytic forms of discrete and continuum hydrogenic wave functions: 


(100|h) 

{nlmi\r) 




TT 


( 3 ) 




{21 + 1)! y 2n(n — / — 1)! \ n / 


n 


-(n-/-l),2/ + 2. 


2Zf 


n 


{Pr\r) =e2pdr ( 1 - — ^ e l,z(p^r + p; ■ r) ) , 

Pr J \Pr 


( 4 ) 

( 5 ) 


are given in atomic units [so that barred quantities f = rruea and pr = Pr'/(uie'u)], where 
r(; 2 ) and iEi{a,h,z) are the Gamma and confluent hypergeometric functions, respectively. 
By the Nordsieck integration techniques jl6Hl9], we can calculate matrix elements of the 
transition operator where if denotes the three momentum transfer, analytically. The 
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response function relevant for transitions to discrete states is found to be 




mi 




= (2/ + l)lli, 

_l)n-z-i j r(2/ + 2) r(2/ + 3) K ; 


( 6 ) 


4n2(2/ + l)\\j {n-I- 1)! r(n + / + 1) T{1 + 3/2) MZ' 

n—l—1 r — —l—2~ 


d 

^ ' Jt 


(7) 


t=0 


which is dimensionless. The response function relevant for transitions to continuum is 

-O 

z,ivi ^' ■ 
r nr/ 

2®ZV(3/t^+P^ + ^^)exp 


fl(”">(K) = j d%\ (p-|e«-'|ls> , 

2Z f 2Zpr 


- 2 ' 7 rZ 

3mea2((K + p^)2 + Z2)3 ((k — Pr)2 + Z2)3(l — exp Pr ) 


Pr=^HT-B-&) 


( 8 ) 


where the factor of mecr^ in the denominator gives the dimension as the energy conserving 
delta function is included in the definition. 

Using the generic response functions obtained above, the single differential cross section 
with respect to energy transfer, da/dT, can be compactly expressed. For elastic scattering 
or discrete excitation to the final discrete level {nl) 


dT 


- 7 mn 

,(e) 27r vl 

with (f - 



2Mh(T - [Eni - Eu)). 


(9) 

( 10 ) 


The factor in Eq. (j^ comes from two sources with each one contributing l/v^\ (i) 
division by flux in da/dT and (ii) the integration of DM scattering angle cosd with respect 
to the energy conserving delta function. Note that it does not lead to physical sigularity 
when taking an extremely NR limit —)■ 0, since the DM flux and kinetic energy both 
approach zero. The magnitude of q is determined by energy conservation, or equivalently, 
the scattering angle cos 6 is fixed once the energy transfer T for such 2-to-2-body scattering 
is known. For ionization processes: 
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d^iion) 


dT 




1 rUy, 

= - -k2 

2tt Uv 


dcos9 


,(e) 


(k= -^q), 


Trip 


with min < 1 , max 


■ 1 , 


kl + kl- 2M{T - B) 
2kik2 


( 11 ) 

< cos 6 * < 1 , ( 12 ) 


where ki = rriyVy and k 2 = {rri^v'^ — are the magnitudes of the initial and hnal 

momentum, respectively. Because the hnal atomic state has two bodies to share the trans¬ 
ferred energy and momentum, the DM scattering angle cos 9 now can span a hnite range for 
a given energy transfer T. 

Next we consider the d^^'* term. Because its Lagrangian differs from term only by a 
kinematic factor 1/q^ (which only cause a rescaling of transition matrix elements), it can 
easily be calculated by replacing to d^f*/q^ in both Eqs. (9) and (11). If both terms 
exist, then one has to take their coherent interference into account, so that the coupling in 
front of the response functions becomes -f d^i^/q"^. 

Unlike the and d^^'^ terms, which are independent of the spins of the DM and the 
scattered particles, the and d!f^ terms give rise to what typically called spin-dependent 
interactions. Their matrix elements for unpolarized scattering involve additionally initial 
spin average and hnal spin sum. For a spinor with spin quantum number s and rris, it yields 


EE 

m' 


{s. m.ISJs. m,) = -s(s + l)4j 


(13) 


With the DM spin Sy and the electron spin Sg = 1/2, the spin averages and sums applied 
to both the DM and electron parts yield a product: Sy{sy -f l)/4. Other than this factor, 
the rest of spatial matrix elements are completely the same as in the and d^i'^ case. As 


a result, the corresponding scattering formula can be obtained by changing 


+ df/q^ 


to |s^(s^ + 1) -f /q^ in both Eqs. (9) and (11). 

It is worthwhile to point out here that there is no interference between the spin- 
independent interactions with ci, di and the spin-dependent one with C 4 , d^, in unpolarized 
scattering, since the trace of a spin matrix is zero. 

Now we consider the cases when DM scatters off the proton instead of the electron. 
Besides the trivial change of LECs, the most important difference is due to the fact that the 
proton is much closer to the atomic center of mass than the electron. After factoring out the 
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center-of-mass motion, the resulting atomic transition operators in its intrinsic frame are 


= e'(14) 
p(p)(^ = (15) 

for the electron and the proton, respectively. This leads to the following change of the 
corresponding response functions: 

Mon) ^ _ Ag) ( 16 ) 

rUe 

Mion) ^ Jl{nl,ion)r^ ^ Ag) ^ ( 17 ) 

^ rUp 

and similarly in the differential cross section formulae, Eqs. 0 and ( pT| . 

Finally we can summarize the above derivation and obtain the differential cross section 
formulae for DM scattering off the hydrogen atom at LO. For transitions to discrete states: 



and for ionizations: 


d^iion) 

dT 


LO 



+ + 1 ) 


M) 

Jf) , ^4 
^4 0 

1 

qz 



p2=2/i(T-B-g2/(2M)) 


R 


(ion) 

f 


(19) 


In these formulae, there are two new response functions defined, which describe the 

nontrivial interference between the spin-independent y-e and y-p amplitude; they are 


p(nO 

^ep 




( 20 ) 


and 


= / d^p^ {pr\R^^^'"\ls) {Pr\e 6{T - B 


-Q. 

q Pr 


2M 2p' 


( 21 ) 























Note that there is no interference between the and terms, nor between the and 
terms, as they involve different spin operators with each of them traceless. 

Even though we mainly concentrate on the LO interaction with DM in this article, we 
shall also consider a few terms at next-to-leading order (NLO); 


42^°^ = ■ <?/) + 4^• ?X)(4/) 

f=e,p,n 

+d[o ■ Qf) + ■ gx)(4/)| + ■ ■ ■ (22) 


They translate into the ones of |13] by cn —)■ h^, cio —)■ hg, dn —)■ li, and dio —)• l 2 - Because 
the spin operators that come with the and terms are mutually independent (after 
spin average and sum) with each other except for the interference between and and 


also with all LO terms, including them into Eqs. (18) and (19) is straightforward by 



where similar argument is applied to the d 


(e,p) 
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and 


terms. 


III. RESULTS AND DISCUSSIONS 

In this section, we give numerical results for two different DM masses: (i) = 1 GeV 

and (ii) = 50 MeV, and with a nonrelativistic velocity = 10“^. Our main purpose is to 
illustrate and discuss how electron and proton respectively contribute to the scattering pro¬ 
cesses for specified EFT interaction terms and reaction channels. For clarity in presentation, 
we ignore all interference terms and assume the DM interaction strengths with electron and 
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proton are the same when making comparisons between the electronic and nuclear contri¬ 
butions. It should be borne in mind that the total cross section is a sum of all contributions 
from electron and proton with interference terms included. 


A. LO Interactions of ci, di, C 4 , and di 

The upper two panels of Fig. [^show the differential cross sections da/dT for DM scatter¬ 
ing with = 1 GeV and the Ci-type interactions. For the nuclear part (right panel), the 
elastic scattering dominates all other channels by orders of magnitude. The reason is obvious: 
Since the momentum scale that determines the nuclear response Kp = ~ ~ 

0.5 keV is smaller than the inverse of atomic size, meO; ~ 3keV, it is a good approxima¬ 
tion that the nuclear charge operator can be expanded as p{Kp) ~ 1 — iHp ■ r + ... Unlike 
the elastic scattering, all inelastic channels have no leading-order contributions because of 
wave function orthogonality. Also because the next-leading-order operator is a dipole oper¬ 
ator, the excitations to final p-orbitals or continuum are more probable than other discrete 
states. For the electronic part (left panel), the results change dramatically. First, as the 
momentum scale that determines the electronic response ~ ~ 1 MeV 

is much bigger than the inverse of atomic size, the electric charge operator p{Ke) = 
becomes highly oscillating. As a result, the elastic differential cross section shows the fa¬ 
miliar form factor suppression. On the other hand, in discrete excitations, one does observe 
much larger cross sections in near-threshold regions. The reason is most of the energy 
transferred T by DM is given to internal excitation; this leaves the 3-momentum transfer 
q = — {Eni — Els)) becoming quite small so that the form factor suppression is 

less severe. Among all reaction channels arising from the c^'^^-type interaction, the ionization 
channel is the dominant one in most of the kinematic region, for it can access more of the 
kinematic phase space with small q. In addition, the peaks near discrete excitation thresh¬ 
olds also provide good observation windows for not only the large cross section but also the 
clean signal of deexcitation photons of definite energies. 

In the lower two panels of Fig. the results for the di-type interactions are shown. 
The numerical changes from the previous Ci-type results are mainly due to adding 1/g^ 
factor in the scattering amplitude (or 1/q^ in the double differential cross section). Because 

= 2mi{T in elastic scattering, this leads to an extra 1/T^ dependence in da/dT and 
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10 ^ 


cf’=1/MeV^; m^=1GeV 


cj'’>=1/MeV^ m^=1GeV 






FIG. 1. Different channels of DM scattering with m-^ = 1 GeV and interaction of = 1/MeV^ 
(upper-left), = 1/MeV^ (upper-right), = 1 (lower-left), = 1 (lower-right). 


can be best illustrated by observing the difference of the (almost) flat line and the power- 
law-decreasing line in the and plots respectively. For discrete excitation channels, 
except for near-threshold region, one expects similar 1/T^ dependence when T gets bigger 
than excitation energies. The case of ionization channel is more intricate, as is to be 
integrated over a range allowed by kinematics; there is no simple scaling from the Ci- to di- 
type results. Overall, the long-range interaction yields sharper energy dependence of da/dT 
than the contact one for all scattering channels considered. The elastic scattering is still the 
best channel to constrain d!f \ and discrete excitations at thresholds and ionization the best 
for d^f’. 

In Fig. 1^ similar plots but with = 50 MeV are shown. The most noticeable differences 
from Fig. [^are (i) the NR DM kinetic energy is smaller so Tmax = = 25 eV and (ii) 

in elastic scattering and discrete excitations, the maximum allowed energy transfers are cut 
off at smaller values: 4.8, 14.0, and 15.7 eV for hnal n = 1, 2, 3, respectively. The latter is 
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due to maintaining energy and momentum conservation in the final two-body system (the 


DM particle and the atom) with < Mjj. Using Eq. (10) and setting the maximum DM 


scattering angle cos 6 = — 1, one can get an approximate formula 


rpinl) _ rji I 

cut / , 71 /T \ o max “T 


AtTI^Mh rn . 


Mh + m. 


{Enl — Els) , 


(24) 


which yields the correct cut off energies just pointed out. Except the cutoffs in energy trans¬ 
fer, da/dT^s of elastic scattering and discrete excitations are the same for both = 50 MeV 
and = 1 GeV in the allowed range of T. Because the associated response functions, which 
depend on the 3-momentum transfer q, are fixed by T and excitation energy, the indepen¬ 
dence of is thus understood. Note that these cutoff energies limit the ability of direct 
LDM searches because the recoil energies are too small to be detected. 


c{°>=1/MeV^ /n^=50MeV 



1; m^=50MeV 



FIG. 2. Different channels of DM scattering with i 
(upper-left), = 1/MeV^ (upper-right), = 

blue curves are the ionization results with = 
Fig. 0 these are linear-log plots.) 


cl'’>=1/MeV^ m^=50MeV 




= 50 MeV and interaction of = 1/MeV^ 
. (lower-left), = 1 (lower-right). The thin 
1 GeV, shown for comparison. (Note. Unlike 
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On the other hand, da/dT of the ionization channel shows different features. First, as 
there are three bodies in the final states (the DM particle, the ionized atom, and the ejected 
electron), energy and momentum conservation does not introduce a kinematic cutoff so T can 
extend to the end point energy Tmax- For this reason that the ionization channel should be 
considered as the golden mode to LDM direct searches. Second, the value of q does depend 


on via Eq. (12); as a result, da/dT is not independent. To make the comparison 
clear, the results of = 1 GeV are plotted with thin solid curves in the same figures. 

At T 17 eV, one observes discontinuities in da/dT on the right panels of Fig. This 
is a combined result of two ingredients: (i) The scattering angle is bounded by Eq. ([T^. At 
T ^ 17 eV, the maximum scattering angle 180° is reached (this energy can also be anticipated 


from Eq. (24) with the excitation energy Eni — Eu being replaced by the binding energy 
—Els ), so the integration range in Eq. ceases to increase for T > 17 eV. (ii) The 
nuclear response function is bigger at backward angels than forward angels, so the 

integral Eq. ( [IT| sensitively depends on the integration range and its discontinuity. On the 
contrary, the electronic response function Re°^'^ is only significant at small angels, therefore 
the discontinuity in the integration range of Eq. (0 does not yield observable results on 
the left panels of Fig. 

The hydrogen atom only has one electron and one nucleon, so the contributions from the 
C 4 - and (i 4 -type interactions are related to the ones of the Ci- and di-type interactions simply 
by a rescaling of one-body spin matrix element as discussed in the last section. We shall not 
repeat these plots, but just note that for other atoms with more electrons and nucleons, the 
spin-dependent cross sections from interaction terms like C 4 and d^ do not receive many-body 
enhancement compared with the spin-independent interactions terms like ci and di. 

As DM interactions with electrons and nucleons are both included in our calculations, 
it is interesting to compare their contributions. Assuming the same coupling constants, 
c^//^ = = 1/MeV^, d^//^ = d^'^ = 1, the comparison shown in Fig. |^gives several important 

features: 


1. In elastic scattering, the nuclear contribution dominates, and is bigger than the elec¬ 
tronic part by several orders of magnitude. Therefore, elastic scattering is not likely 
to be a good channel of constraining the LO DM-electron interactions, if the LO 
DM-nucleon interactions are present and not unnaturally suppressed. 
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FIG. 3. Comparison of DM cross sections with the electron (xe) and proton (xp) iii ^ hydrogen 
atom for selected channels including (i) elastic (Is), (ii) discrete excitation to 2s, and (iii) ionization 
(ion). The interactions are taken to be (Left) = 1/MeV^ and (Right) = 1. 

Interference terms due to ye and XP amplitudes are ignored. 


2. In discrete excitations, the nuclear and electronic contributions have sharp crossovers 
at energies slightly bigger than excitation energies. If a detector is able to resolve 
these peaks where electronic contributions clearly dominate, then it can be useful for 
setting more stringent limits on the LO DM-electron interactions. 

3. In ionization processes, unlike discrete excitations, the electronic contributions gener¬ 
ally dominate over the nuclear parts up to some T beyond the ionization thresholds. 
As a result, the LO DM-electron interactions can hopefully be constrained in broader 
kinematic regions. 

In Fig. 1^ we study the m^-dependence of the crossover energy below which the DM-electron 
cross section begins to be bigger than the DM-proton one (assuming the same coupling 
strength) via the Ci- or di-type interaction that gives rise to hydrogen ionization. The hrst 
thing to notice is in both types of interactions and the considered range of (50 MeV to 5 
GeV), there exist certain ranges of DM energy transfer T where the ionization processes are 
more sensitive to the DM-electron interaction than the DM-nucleon one. Furthermore, one 
observes that the crossover energy for the di-type interaction is larger than the one for the 
ci-type interaction. The main reason is the 1/g^ factor appearing in the double differential 
cross section gives more weight to the response function at low which enhances the role 
of electrons on one hand and suppresses the role of nucleons on the other. 
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(MeV) m,i (MeV) 


FIG. 4. Energy transfer of DM, T^s, below which scattering with electron yield bigger da/dT 
than proton in the hydrogen ionization (assuming the same xe and XP coupling strengths), plotted 
against DM mass (Left) for the ci-type and (Right) for the di-type interaction terms. 

Therefore, it is reasonable to conclude that the best observational window to look for 
the LO DM-electron interactions is ionization processes near threshold, in particular for 
LDM with < Mh- The discrete excitation peaks (which need good energy resolution 
of detectors) also provide good supplements. Although hydrogen can hardly be a good 
candidate for detecting LO DM-electron interactions for the low energy transfer T ~ 10 — 
20 eV is far below the current detector thresholds, however, our above conclusion makes good 
sense for practical detector species made of heavy atoms: Not only the ionization thresholds 
of atomic inner orbitals can be as high as a few or tens of keV which are observable in current 
detectors, but also there are more than one ionization peaks which can provide additional 
information. 

B. NLO Interactions of cn, dn, cio, and dio 

Consider now the interaction terms of cn and dn, the results for = 1 GeV are pre¬ 
sented in Figs. The main change from the corresponding plots of Ci and di is the extra 
factor appearing in the differential cross sections. For elastic scattering or discrete excita¬ 
tions away from threshold, this factor introduces an extra dependence on T, and again can 
be best seen by a comparison of the and curves for elastic scattering. For ionization, 
the impact of the factor, which is to be integrated over some allowed range, however can 
not be easily factored out. 
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FIG. 5. Different channels of DM scattering with m-^ = 1 GeV and interaction of = 1/MeV^ 
(upper-left), = 1/MeV^ (upper-right), = 1/MeV (lower-left), = 1/MeV (lower-right). 

Our previous argument that da/dT^s of elastic scattering and discrete excitations are 
independent of also applies to the cases of cn and du (also cio and dio to be discussed 
later). Therefore, we do not repeat the results for = 50MeV and just point out that 
they are the same as for the = 1 GeV case in the regions bounded by the cutoff energies 


given by Eq. (24). Also, while the ionization processes do have ruy dependence, it does not 


differ from what has been shown in Fig. [^for the case of Ci and di in a significant way. 

Similarly to the Ci- and di-type interactions with DM, elastic scattering is the best to 
constrain the and d^^^ terms, while inelastic scattering at discrete excitation peaks and 
of ionization are more suitable for the and d!fl terms. To further disentangle the de¬ 
pendence of da/dT on ci, di, cn, and du, the scaling of da/dT with energy transfer T can 
provide useful guidance: For example, in elastic scattering and discrete excitations away 
from thresholds, the energy dependence of da/dT on the |cip, |dip, |ciip, and |diip terms 
is T°, T“^, T^, and T~^, respectively. 
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The pattern regarding the competition of electronic and nuclear contributions in discrete 
excitations and ionization with the Cn and dn terms is similar to the case with the Ci and 
di terms; sharp crossovers in discrete excitations and some ranges of electronic dominance 
near ionization threshold. Some examples are given in Fig. 



FIG. 6. Comparison of DM cross sections with the electron (xe) and proton (xp) iii ^ hydrogen atom 
for selected channels including (i) elastic (Is), (ii) discrete excitation to 2s, and (hi) ionization (ion). 
The interactions are taken to be (Left) = 1/MeV^ and (Right) = 1/MeV. 

Interference terms due to ye and XP amplitudes are ignored. 


The energy transfer Txs below which the electronic contribution is bigger than the nuclear 
one, assuming = 1/MeV^ and d^^l = d^fi = 1/MeV, in the ionization processes 

is plotted in Fig. against Notice that the values of for the Cn- and dn-type 

interactions at a given are both reduced in comparison with the cases of the LO Ci and 
di terms, respectively. This is in agreement with the expectation that the extra factor 
in the double differential cross section reduces the weight of the small region, so the 
electronic contribution is relatively suppressed than the nuclear part. 

Similar to the LO case, the contributions from the target-spin-dependent cio and dio 
terms can be obtained from the target-spin-independent results of cn and dn, simply by 


adding factors due to spin matrix elements (see Eq. 13). Therefore, all observations and 
conclusions made in the Cn and dn cases apply to the Cio and dio. 

However, regarding the competition between the electronic and nuclear contributions in 
scattering processes involving the cio or dio term, there is a subtlety arising from the natural 
scales of and d^l^ /d[g^. If one takes the point of view that both EFT interaction terms 

of electrons and nucleons are matched to a more fundamental theory at some high scale A, 
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FIG. 7. Energy transfer of DM, Txs, below which scattering with electron yield bigger da/dT 
than proton in the hydrogen ionization (assuming the same xe and XP coupling strengths), plotted 
against DM mass (Left) for the cn-type and (Right) for the dn-type interaction terms. 

then it is reasonable to anticipate the possibility that ~ ^ ~ 

the other hand, the masses of an electron and a nucleon differ by three orders of magnitude. 

If the cio and dio terms are matched to a relativistic theory, for example, (xx)(/* 75 /) and 
some high scale, the resulting nonrelativistic EFT expansion at NLO will 
involve the expansion parameter q/rrif to hrst order. In such cases, then one should expect 
4oVcio\ ~ ~ 2 X 10^. This in turn would largely increase the sensitivity 

of discrete excitation peaks and ionization processes on the NLO DM-electron interaction 
terms such as Cio and dio- 

An example is given in Fig. For < 160 MeV and < 240 MeV respectively for the 
Cio and dio terms, the electronic contributions are larger than the nuclear ones in the entire 
allowed ranges of T < For heavier the crossovers both happen at energies 

further away from ionization thresholds, ~ 50 eV and 100 eV respectively for the cio and dio 
terms - much bigger than other interactions terms previously discussed. 

Briefly concluding this subsection, we point out that the best observational window to 
look for the NLO DM-electron interactions terms including cn, dn, Cio, and dio is still 
the ionization processes near threshold and the discrete excitation peaks. The different 
energy dependence of da/dT from the LO terms in principle provides a way to disentangle 
them. Furthermore, because of the huge difference between the masses of an electron and 
a nucleon, interaction terms that depend on the relativity of scattered particles can be 
further separated. In most situations, such NLO DM-electron interactions can be sensitively 



14.45 
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FIG. 8. Comparison of electronic (solid lines) and nuclear (dotted lines) contributions to da/dT for 
selected DM mass with (Upper) = 1/MeV^ and (Lower) d^^Q = d^^Qm^/ml = 

1/MeV. 

constrained without much background arising from similar DM-nucleon interactions because 
atomic electron can be very relativistic while atomic nuclei and nucleons inside are mostly 
nonrelativistic. 


IV. SUMMARY 

In this paper, we study the scattering processes of sub-GeV DM particles and hydrogen 
atoms, including elastic, atomic discrete excitation, and atomic ionization channels. The in¬ 
teractions of DM with electrons and nucleons are both included and formulated in a general 
framework based on nonrelativistic effective held theory. In addition to the leading-order 
spin-independent and spin-dependent contact terms, we also include the possibility of long- 
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ranged DM interactions and a few next-to-leading-order terms. Some of the interaction terms 
yield orthogonal scattering amplitndes, bnt there are also interference terms. Disentangle¬ 
ment of varions interaction terms can in principle be done by their different dependence on 
DM energy deposition in scattering cross sections. 

On the assnmption of same dark matter conpling strengths, it is fonnd that DM-electron 
interactions dominate the inelastic transitions to discrete excited states and ionization con- 
tinnnm aronnd their threshold regions (sizes of these regions depend on interaction types), 
and DM-nncleon interactions become more important with increasing energy and dominate 
in elastic scattering. These conclnsions can be nsed to gnide the searches of snb-GeV DM 
interactions in optimal experimental conhgnrations and kinematics. For DM-electron inter¬ 
actions, the inelastic peaks of discrete excitations and ionizations in scattering cross sections, 
which can be taken as smoking-gnn signals of DM scattering, can fnrther increase an exper¬ 
iment’s constraining power. For DM-nncleon interactions, althongh the elastic scattering is 
the best channel, however, for light DM particles which can not deposit observable energies 
in detectors, one has to rely on the high energy part of ionization processes. 

The energy and momentnm transfers involved in snb-GeV DM scattering overlap typ¬ 
ical atomic scales, so stndies of issnes snch as binding effects and electron/nnclear recoil 
mechanism, which play important roles in interpreting experimental data, reqnire detailed 
many-body calcnlations. This case stndy of hydrogen, where both binding and recoil can 
be taken into acconnt most simply, therefore provides nsefnl qnalitative nnderstanding of 
what to be anticipated in snb-GeV DM scattering off practical detector materials snch as 
germaninm and xenon. 
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